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We calculate the contribution of superconducting fluctuations above the critical temperature T c to 
the transverse thermoelectric response a xy , the quantity central to the analysis of the Nernst effect. 
The calculation is carried out within the microscopic picture of BCS, and to linear order in magnetic 
field. We find that as T — * T c , the dominant contribution to a xy arises from the Aslamazov-Larkin 
diagrams, and is equal to the result previously o btained from a stochastic time-dependent Ginzburg- 



Landau equation [Ussishkin, Sondhi, and Huse, axXiv:cond-mat/0204484|. We present an argument 



which establishes this correspondence for the heat current. Other microscopic contributions, which 
generalize the Maki- Thompson and density of states terms for the conductivity, are less divergent 
as T -> T c . 



I. INTRODUCTION AND DISCUSSION OF 
RESULTS 

In a superconductor, fluctuations of the superconduct- 
ing order parameter above the transition temperature 
T c affect various properties such as the magnetic sus- 
ceptibility and transport coefficients. The study of su- 
perconducting fluctuations has a long history (for re- 
views, see, e.g., Refs. More recently, interest in 
fluctuation phenomena was renewed with the discovery 
of high-temperature superconductors, where their short 
coherence lengths, high critical temperatures, and lay- 
ered structuresj-imply a large regime for the observation 
of fluctuations n 

One experiment that has aroused particular interest 
recently is that of the Nernst effect: A temperature gra- 
dient (— VT) || x is applied in the presence of a magnetic 
field B || z, and the electric field response (in the ab- 
sence of transport electric current) is measured in the y 
direction. Below T c in the vortex state the Nernst effect is 
large due to vortex motion, while in the normal state it is 
typically very small. In experiments in low temperature 
superconductors no sign of superconducting fluctuations 
was reported as the temperature was raised above T c .a 
In contrast, several different experiments did observe the 
appearance of a fluctuation tail above the critical tem- 
perature in the. fttepfist signal of the high-temperature su- 
percanductoraa-ElDEi (and also in the related Ettinghausen 
effecta). More recently, the Nernst effect above T c has at- 
tracted considerable attention with measurements show- 
ing a sizeable Nernst signal, suell above T c , in particular 
in the underdoped regime.ESEJ 

While a Nernst experiment is carried out under open 
circuit conditions, the transport coefficients which arise 
naturally in a theoretical description are those which re- 
late the transport electric and heat currents to the elec- 
tric field and temperature gradient, 
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Here, a is the conductivity tensor, k a tensor of thermal 
conductivity, and a, a the thermoelectric tensors (which 



obey the Onsager relations, a = Ta). Applying the open 
circuit condition to Eq. ([!]), the Nernst coefficient is ex- 
pressed in terms of the conductivity and thermoelectric 
tensors, 
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The transverse thermoelectric response a xy , the quantity 
on which this paper is focused, is of primary interest for 
understanding the effect of superconducting fluctuations 
on the Nernst signal (as discussed below). 

In a recent paper, Ussishkin, Sondhi, and Huse dis- 
cussed the contribution of superconducting fluctuations 
to the thermoelectric and thermal conductivity ten- 
sors using a stochastic time-dependent Ginzburg-Landau. 
equation (TDGL) in the limit of Gaussian fluctuations.!! 2 ] 
In this paper we revisit the calculation of the transverse 
thermoelectric response a xy using a diagrammatic calcu- 
lation within BCS theory. The details of this calculation 
are presented in subsequent sections. In the remainder 
of this section we present and discuss the results of this 
paper. 

We calculate a xy above the critical temperature T c , to 
linear order in the magnetic field B || z, and to leading 
order in T — T c . We find that, in two and three dimen- 
sions, the contribution of superconducting fluctuations to 
the transverse thermoelectric response is 
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is the magnetic length, and 
£(T) oc (T — T c )~ 1//2 is the coherence length of the super- 
conducting order parameter. 

It is well known (see, e.g., Refs. [l],^) that super- 
conducting fluctuations enhance the conductivity above 
T c due to rbfjttJji the Aslamazov-LarkinEj and the Maki- 
ThompsonEaia contributions (there are also density of 
states terms, which are less important for the conductiv- 
ity). A similar identification of the microscopic contribu- 
tions applies to other transport coefficients. In the case of 
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the transverse thermoelectric response, as the superscript 
in Eq. (||) suggests, we find that the leading order contri- 
bution to a xy is due to the Aslamazov-Larkin diagrams 
alone. The contribution of the Maki-Thompson and den- 
sity of states diagrams is less divergent as T — ► T c . 

Physically, the Aslamazov-Larkin diagrams correspond 
to the contribution of thermal fluctuations of the order 
parameter. Their contribution to a xy may be viewed 
either as the transport heat current carried by such fluc- 
tuations when they respond to an electric field, or as 
the transport electric current carried by the fluctuations 
as they respond to a temperature gradient (all in the 
presence of the magnetic field). The same physics is 
identically described by the Gaussian approximation to 
a stochastic TDGL. (We will have more to remark on 
the correspondence between the two approaches in subse- 
quent sections.) Indeed, the result obtained in this paper, 
Eq. (||) , is identical to the result obtained in the Gaussian 
approximation to the stochastic TDGL in Rcf. [l2| 

Before discussing the experimental situation, we com- 
ment on two assumptions made in this paper. First, 
we assume that the order parameter has s-wave sym- 
metry. In the context of the high-temperature super- 
conductors it is of interest to consider also the case of 
d-wave symmetry in this approach. We note here that 
this will not affect the conclusions of this paper: the re- 
sults of the stochastic TDGL would still correspond to 
the Aslamazov-Larkin contribution; and the arguments 
showing that Maki-Thompson and density of states terms, 
are less divergent remain valid in this case as well.t£l 
(We do not consider here the related issue, of whether 
nodal quasiparticles, which appear when the tempera- 
ture is lowered and the condensate is formed, contribute 
to ot X y) Second, we assume particle-hole symmetry (i.e., 
neglecting any contributions which arise due to asymme- 
try around the Fermi surface in properties such as the 
density of states). For quantities that do not vanish in 
this limit this is a very good approximation for a BCS 
superconductor. (On the other hand, particle-hole sym- 



metry implies that a xy = 



= 0, and therefore 



in calculating these transport coefficients it is necessary 
to break this symmetry.) We note that the contribution 
of the normal metallic state to a xy also vanishes in this 
limit. However, this is not required by symmetry, and 
indeed is no longer the case once superconductivity is 
taken into account. 

We now return to the discussion of the Nernst coeffi- 
cient v N . As noted in Eq. (g), v N is related to both the 
conductivity and thermoelectric tensors. However, the 
main effect of superconducting fluctuations on the Nernst 
signal above T c is due to a xy . Indeed, the contribution 
of fluctuations to the second term in the numerator of 
Eq. (||) is small due to considerations of particle-hole 
symmetry. Moreover, not too close to T c the conduc- 
tivity is dominated by the normal state contribution. It 
follows that the main contribution of superconducting 
fluctuations to the Nernst signal (to linear order in B) is 
a xy /o~ xx , with a xx being the normal state contribution. 



Since the result for a xy , Eq. (||), depends only on the 
coherence length £(T), and in a simple manner, com- 
parison with experiment should apparently be straight- 
forward. In Ref. such a comparison for a high- 
temperature superconductor was presented. On the 
other hand, in low-temperature superconductors, for 
which BCS theory is certainly applicable, the appearance 
of the fluctuation tail in the Nernst signal was not pre- 
viously reported to the best of our knowledge.!!! The rea- 
son for this is that low-temperature superconductors are 
typically also good conductors in the normal state. Con- 
sequently a xy I a xx , the contribution of superconducting 
fluctuations to the Nernst signal, is strongly suppressed 
in bulk low-temperature superconductors. 

The situation can be improved considerably by look- 
ing at a thin film, which is effectively a two dimensional 
superconductor if the coherence length is larger than the 
film thickness. First, the fluctuation tail of a xy is en- 
hanced by going to lower dimensionality, as is evident in 
Eq. (H|). [The result for two dimensions in Eq. (||) is to 
be divided by the film thickness to obtain the result for 
a thin film.] Second, such films may have a significantly 
lower normal state conductivity. Taken together, these 
effects may considerably enhance the contribution of fluc- 
tuations to the Nernst signal. A similar situation may 
occur in a layered structure with weak coupling between 
the layers and with a low normal-state conductivity, r-as 
is the case for the high-temperature superconductors.113 

In the remainder of the paper we present the details 
of our calculation. In Sec. [il] we discuss the Kubo for- 
mula for a xy , and discuss an important aspect of the 
problem, namely the role of bulk magnetization currents 
and the proper subtraction of their contribution.EJ We 
also briefly present known results for the propagator of 
superconducting fluctuations. The diagrams that appear 
in the calculation of a xv , and the physics they describe, 



are discussed in Sec. [IL In Sec. IV we discuss the cal- 



culation of the heat current vertex. A general argument 
regarding its calculation is given, establishing the corre- 
spondence to the heat current in the TDGL. In Sec. [V| 
we calculate the Aslamazov-Larkin diagrams, and obtain 
Eq. (||). The Maki-Thompson and density of states dia- 
grams are considered in S ec. VI. Finally, we summarize 
our discussion in Sec. 



VII 



II. FORMALISM 



In this section we present a few results which will form 
the basis of the calculations in the sections that follow. In 



Sec. II A we discuss linear response theory for the trans- 
verse thermoelectric response. We present the Kubo for- 
mula for a xy , and discuss the subtraction of bulk mag- 
netization currents. Presented here only for the sake of 



completeness, Sec. [I B briefly discusses known results for 



the propagator of superconducting fluctuations. 
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A. Kubo formula for a xy 

In this paper the thermoelectric tensor is considered 
by calculating the heat current response to an electric 
field. Alternatively, one could consider within linear re- 
sponse theory the electric current response to Luttinger's 
"gravitational field" Ea The result for a xy is of course in- 
dependent of which formulation is used, and results are 
presented in terms of one of them only for convenience. 

The heat current response to an electric field is re- 
lated to the heat current-electric current correlator by 
using the standard Kubo formula. Here, we are inter- 
ested in the calculation of this response to linear order in 
the magnetic field. This amounts to introducing an acL 
ditional current vertex coupled to the magnetic field.t£l 
Here we present this result as the Kubo formula for the 
linear response to both an electric and magnetic field. 
The electric field E || x and the magnetic field B || z are 
introduced at finite frequency and wavevector, respec- 
tively, using the vector potential 
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The heat current in the y direction, in response to the 
electric and magnetic fields (in the d.c. limit), is given by 

^ = -n!^oi Re [ A ^ n ^^- +i o. (5) 
Here, the three current correlator A is defined by 
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where f2 m = 2nmT is a bosonic Matsubara frequency 
(with units in which h = k B = 1), and the upper limit of 
integration over imaginary times r and t' is the inverse 
temperature (3 = 1/T. In Eq. (||), an analytic continu- 
ation of Q m to real frequencies is performed before the 
zero frequency limit is taken. 

An important aspect of the calculation of the trans- 
verse thermoelectric response,., discussed in detail by 
Cooper, Halperin, and RuzinJiH is the need to account 
for bulk magnetization currents. This issue arises because 
the microscopic electric and heat currents, as calculated 
by the Kubo formula, are composed of transport and 
magnetization currents, 
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The magnetization currents are currents that circulate 
in the sample and do not contribute to the net currents 
which are measured in a transport experiment. On the 
other hand, they do contribute to the total microscopic 
currents, and it is thus necessary to subtract them from 
the total currents to obtain the transport current re- 
sponse. In the presence of an applied electric field, it 



was shown in Ref. [l?] that the magnetization current is 
given by 



Jmag = CM X E, 



(8) 



where M is the equilibrium magnetization (in the absence 
of the electric field) . It then follows that the transverse 
thermoelectric response is given by 
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hi 



(9) 



where jy/E is found using the Kubo formula, Eqs. (||) 
and (§). 

In calculating the electric current response to a "gravi- 
tational field" ip, a similar situation arises, where in order 
to obtain the transport current the electric magnetization 
current-has to be subtracted. The latter is given in this 
case bytll 



= -cM x VV>. 



(10) 



In Ref. |l2| the total heat current response to an electric 
field and the total current response to a temperature gra- 
dient were calculated using the stochastic TDGL. The ap- 
parent discrepancy between the results for the total cur- 
rents and the Onsager relations was invoked to demon- 
strate the need of subtracting out the magnetization cur- 
rents. In contrast, in the linear response formalism the 
calculation for the electric current yields its response to a 
"gravitational field" gradient, and this apparent discrep- 
ancy does not arise. The magnetization currents (and 
the total currents) trivially obey the same Onsager re- 
lations obeyed by the transport currents. Nevertheless, 
the magnetization currents must be subtracted to obtain 
the correct result for the transport coefficients. 



B. Fluctuation propagator 

The contribution of superconducting fluctuations to 
the current correlator in Eq. (^) is calculated in this pa- 
per for a BCS superconductor (with s-wave symmetry), 
with Hamiltonian 

k.cr k,q,<7 

+ A C k',T C -k'+qa C -k+q,iCk,T- (11) 

k,k',q 

Here, £k = k 2 /2m is the kinetic energy of the electrons, 
a =T>i is their spin, Uq is the disorder potential (with 
the usual Gaussian distribution), and A < is the attrac- 
tive BCS interaction (where only states with energy dif- 
fering from the Fermi energy by at most lo d participate in 
the interaction term). The relevant diagrams for super- 
conducting fluctuations are calculated using the finite- 
temperature diagrammatic technique. This approach is 
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q, u 




FIG. 1: Diagram for the non-interacting two particle propa- 
gator II(q, cum)- Lines with arrows are electronic Green func- 
tions, and disorder is shown in the ladder approximation by 
dashed lines. 



analogous to the one used, e.g., in the case-.of the con- 
ductivity, leading to the Aslamazov-Larkir£3 and Maki- 
Thompsontala contributions. (For a detailed account 
of the diagrammatic calculation of the conductivity, see, 
e.g., Ref. |[) A basic ingredient in this calculation is L, 
the propagator of superconducting fluctuations. 

Accounting for the electron-electron interaction in the 
ladder approximation, the fluctuation propagator L is 
related to the non-interacting two-particle propagator II 
through 



£(q,w m ) = [A 1 -II(q,cJ m )] 



(12) 



(see Fig. [I] for the diagram of II) . To obtain the retarded 
fluctuation propagator, the Matsubara frequency is ana- 
lytically continued to the real axis (iuj m — * uj + iO) and 
II is calculated to leading order in q and u>. Assuming 
particle-hole symmetry, the retarded fluctuation propa- 
gator is then 



L*(q >W ) 
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(and the advanced fluctuation propagator is L j4 (q, cj) = 
[L R (q, u>)]*). In Eq. (|l3|), v is the density of states per 
spin, e = ln(T/T c ) w (T - T c )/T c is the reduced temper- 
ature, t bcs = 7r/8T c , and 
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Here, r e i is the elastic scattering time, D = v^r^/d is 
the diffusion constant (for d dimensions), and ip(x) is 
the digamma function. The parameters in Eq. (O) are 
directly related to the coefficients appearing in a TDG L 
for the order parameter. In particular, £(T) = \frjj~e is 
the superconducting coherence length, and t bcs is the 
relaxation time for the order parameter fluctuations. 



III. DIAGRAMS AND INTERPRETATION 

In this section we present the diagrams which appear in 
the calculation of the correlator ([}]), and discuss the phys- 
ical processes which they represent. For this purpose, the 
microscopic picture is perhaps best recast in terms of a 



quantum functional integral approach. We begin this 
section by briefly discussing this approach. 

The expectation value of a current operator j (which 
can be either the electric current or the heat current) in 
response to a driving field <p may be expressed in terms 
of an imaginary time functional integral 
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Here, ip and ip are the fermion fields, and the action S is 
given by 



S 







dr / g?x 



^^ip a {x)d T ip a {x) +H(a 



(16) 



where (3 is the inverse temperature, x = (x, r), and TC is 
the Hamiltonian density. Introducing a pairing field A 
via the usual Hubbard-Stratonovich transformation, the 
expectation value of the current may be rewritten as 
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Here, S c s is the effective action for the pairing field which 
is obtained by integrating out the fields ip£2 and 
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where So is the part of the action S which is quadratic 
in ip. The calculation of the contribution of supercon- 
ducting fluctuations to the current proceeds by applying 
a Gaussian approximation to Eq. ( |l7| ) . More specifically, 



and (j) 



AA0 



to 



this involves expanding both SofffA, A 
second order in A and A. 

Consider first the calculation of the conductivity, in 
which case <p is the electric field. In Eq. ( |l7| ) the field ap- 
pears in two places, namely in ( j ) A ^ and in the effective 
action S e s[A,A,<p]. The Aslamazov-Larkin approxima- 
tion involves keeping the field dependence in the effective 
action only. The resulting expression is equivalent to a 
calculation using a stochastic TDGL also done at Gaus- 
sian order: The quantity ( j } AA is the current associated 
with the order parameter configuration. The response of 
the order parameter to the field is described by the effec- 
tive action, which is identical to the TDGL description. 
In particular note that in the TDGL the electric field 
is coupled to linear order to the current associated with 
the order parameter, as is described by the Aslamazov- 
Larkin diagram. Not included in the Aslamazov-Larkin 
approximation are the terms obtained by keeping the 
field in (j)AA0- These describe corrections to the nor- 
mal state response modified by the presence of the order 
parameter, and are the Maki- Thompson and density of 
states corrections. 

The situation for the transverse thermal response is 
somewhat different than it is for the conductivity as we 
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FIG. 2: Diagrams arising in the calculation of the transverse thermoelectric response [see Eqs. (^)-(Q)]. A wavy line represents 
the fluctuation propagator, and the lines with arrows are the electronic Green functions. The vertices are the heat current 
vertex jy (open circle), the electric current vertex coupled to the electric field j% (full circle), and the electric current vertex 
coupled to the magnetic field jt (no circle). "Mirror image" diagrams, which may be obtained by reversing the arrow direction 
on each of the electronic Green functions in the diagrams above [except diagrams (b), (d), and (f), for which the "mirror image" 
is not a new diagram], are not presented [See Fig. H below for the "mirror image" of diagram (a)]. Also not presented are the 
different possibilities of adding disorder to each of these diagrams. 



are considering the linear response to two fields (since are 
considering the heat current response to both electric and 
magnetic fields). The resulting diagrams are presented in 
Fig. H (in most cases, the diagrams have "mirror images" 
which are not presented in the figure, but are of course 
also taken into account). Before proceeding a word on 
nomenclature: while the situation here is a bit different 
than in the case of the conductivity, we will refer to the 
diagrams corresponding to the TDGL contribution as the 
Aslamazov-Larkin diagrams; all other diagrams will be 
collectively referred to as Maki- Thompson and density 
of states diagrams, although they do not correspond to 
corrections to the normal state transverse thermoelectric 
response only, as discussed below. 

The diagram in Fig. |2|(a) and its "mirror image" are 
the Aslamazov-Larkin diagrams, which correspond to the 
contribution of the stochastic TDGL. To obtain these di- 
agrams, the electric and magnetic fields are retained in 
the effective action in Eq. ([Tt]) only, and the average over 
the current operator (which in this case is the average 
over the heat current operator, (J q )aa) gives the heat 
current associated with the order parameter configura- 
tion. Moreover, the motion of the order parameter de- 
scribed by these diagrams is that of the TDGL, with the 
electric and magnetic fields coupled to linear order to the 
electric currents associated with the order parameter con- 
figuration. This correspondence will be revisited below: 
In Sec. IV we discuss the heat current associated with 



we calculate the Aslamazov-Larkin diagrams for a xy and 
find that they give the same contribution as that found 
using a stochastic TDGL in Rcf. [l^. 

The rest of the diagrams describe a variety of processes 
which involve corrections to normal state properties, and 
may be understood along similar lines. Diagrams (b) 
and (c) of Fig. describe a correction to the normal 
state thermoelectric response due to the order parame- 
ter fluctuations, with Maki-Thompson [diagram (b)] and 
density of state [diagram (c)] contributions, but with the 
order parameter responding to linear order to the mag- 
netic field. Likewise, diagrams (d) and (e) describe the 
normal state response to a magnetic field in the presence 
of superconducting fluctuations affected by the electric 
field. Diagrams (f) and (g) describe the heat current 
associated with an order parameter configuration (as in 
the Aslamazov-Larkin diagrams), but with the dynamics 
of the order parameter modified by a term not captured 
by the TDGL. (It is interesting to note that in the mi- 
croscopic picture there are corrections to the TDGL in 
the order parameter motion in this case.) Finally, dia- 
grams (h)-(m) describe corrections to the normal state 
transverse thermoelectric response ot X y 



IV. RELATION BETWEEN CURRENT 
VERTICES 



the motion of the order parameter ( j Q ) AA and its con- In this section we consider the calculation of the trian- 
nection to the heat current in the TDGL, and in Sec. |v| gular block appearing in the Aslamazov-Larkin diagram 
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[see Fig. ||(a)]. In this diagram, the microscopic current 
vertex j can be either an electric current vertex j e or a 
heat current vertex j Q . This diagram corresponds to the 
current in the presence of an order parameter configu- 
ration, and is thus directly related to the current which 
appears in the TDGL. Accordingly, while the microscopic 
current vertex is denoted with j, we denote the current 
vertex presented by the diagram in Fig. [|(a) with J. 

The result for the electric current vertex J^-js well 
known (and is needed, e.g., for the conductivity^) . Our 
main concern here is with the heat current vertex J Q , for 
which we establish the following result: At Q = TO = 
[for conventions regarding incoming and outgoing ener- 
gies and momenta, see Fig. ||(a)], the electric and heat 
current vertices are related by 



q+Q/2, 



2e 



(19) 



Heuristically this form is expected for a preformed pair 
of charge — 2e; but this is the BCS limit, for which an 
explicit calculation is needed. Together with the known 
result for J e [see Eq. (22) below] this allows the calcula- 
tion of the Aslamazov-Larkin diagrams in Sec. [v[ as well 
as obtaining the expression for the heat current in the 
TDGL. £3 

For completeness, we consider first the calculation of 
the electric current vertex J e . The vertex J e is needed 
in the calculation of the Aslamazov-Larkin diagram to 
leading order in external frequencies and momenta, and 
it is thus sufficient to set the external frequencies to zero, 
Q m = Lu m — 0, and consider linear order in the wavevec- 
tors q and Q. There is no linear term in Q: indeed, it 
is straightforwarid-|to show that 3 e (q = £l m = uj m = 0) is 
symmetric in Q.Ea In the calculation of the electric cur- 
rent vertex we thus set Q = fl m = uj m = 0, and calculate 
the vertex to linear order in q, J e (q). 

Setting Q = Q m = u> m = in J e , the diagram for 
the electric current vertex is^as in Fig. |^(b), but with 
uj m = (and j = j e = — eVk).E-3 The same diagram may 
be obtained by inserting an electric current vertex in the 
diagram for II(q, oj m = 0), given in Fig. ^ (and correctly 
accounting for the spin indices). Using the relation 



V k G(k, e m ) = v k G(k, e r< 



(20) 



where G(k, e m ) is the electron Green function, the fol- 
lowing result is obtained, 



J e = -2eV q n(q,w m = 0). 



(21) 



Note the appearance of the Cooper pair charge, — 2e (de- 
noted below as e*). The expansion of Eq. (E1J) to linear 
order, in q gives the result for the electric current ver- 
texJla 



J e (q) = -4erjvq 



(22) 



{rj and v were defined in Sec. |II B| ) . The familiar ex- 
pression for the electric current in the TDGL (in the ab- 
sence of fields), follows from this result. More precisely, 
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FIG. 3: Diagrams used in the text for evaluation of the cur- 
rent vertices J e and J Q appearing in the Aslamazov-Larkin 
diagrams. Lines with arrows are electronic Green functions, 
and dashed lines denote disorder, shown here in the ladder 
approximation. The microscopic current vertex j represents 
either the electric current vertex j e = — evk or the heat cur- 
rent vertex j Q = ie m Vk- 



Eq. (|2^) is the electric current associated with a pairing 
field configuration A = e lq r , using 



V [A*(-iVA) + c.c] 



(23) 



[The conventional TDGL form, with I /2m* replacing r\v 
(as in, e.g., Ref. [l]) is obtained after rescaling the pairing 
field.] 

We now reconsider the calculation of the electric cur- 
rent vertex, this time with arbitrary w m , as a first step 
towards establishing Eq. (p^). Setting Q = Q m = 0, the 
electric current vertex J e is presented in Fig. |^(b) (with 
j = j e = — evk). Using Eq. (p0|), the expansion of J e 
to linear order in q is equivalent to introducing a sec- 
ond velocity vertex. The vertex J e is then written as q 
times the contribution of square diagrams as in Fig. ||(c) 
(with j = j e ). The following observation follows from 
the structure of these diagrams: To linear order in q, 
but at arbitrary u) m , the electric current vertex has the 
structure 



0J„ 



(24) 
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where the function / results from integration over all in- 
ternal momenta in the diagrams. The function / depends 
only on the energy variables appearing in the electronic 
Green functions, and running along the two sides of the 
diagrams (we use here the fact that disorder scattering is 
elastic). In addition, because of the symmetric structure 
of the diagrams in Fig. ||(c), we have 



f(er, 



(25) 



The structure of the result for J e as presented Eqs. (|24|) 
( p5| ) is sufficient for establishing Eq. jl9| ) to linear order 
in q; it is unnecessary to evaluate / explicitly. 

As with the electric current vertex, for the calculation 
of the Aslamazov-Larkin diagram we only need the heat 
current vertex J Q to leading order in wavevectors and 
frequencies. As it turns out [cf. Eq. jl9|)], this is one or- 
der higher than the leading order in J e . By symmetry 
of the structure of the diagrams for the current vertices, 
they are invariant under Q, tt m — > — Q, — Q rn . For the 
expansion of J Q , this shows that there is no term linear 
in Q or fl m alone, but does not exclude a term propor- 
tional to Qf2 m . In the calculation below, we use the fact 
that in the Aslamazov-Larkin diagrams [Fig. ||(a) and its 
"mirror image"] we need the heat current vertex in the 
y direction, perpendicular to Q || x, and thus do not 
consider such a term. 

We thus proceed by setting Q = Q m = 0. The heat 
current .uertex J Q is then given in Fig. |^(b) (withj = j Q = 
*emVk)£-3 As in the case of the electric current vertex 
J e , the expansion of J Q to linear order in q amounts to 
the introduction of a second velocity vertex, resulting in 
square diagrams as in Fig. ^(c) (with j = j Q ). It follows 
from the structure of these diagrams that to linear order 
in q (and at arbitrary oj m ) the heat current vertex has 
the structure 



J Q = qj^ ie m f{e r , 



+ W m ), 



(26) 



with Eq. (pi]), we find the relation between the current 
vertices, Eq. (|l9|). 

We note that to this order of the calculation J Q does 
not have a branch cut after analytic continuation of iu> m 
to the w-plane. After analytic continuation, we thus have 
to linear order in q and u> 



2e 



J e (q) = 2i]i'ujq. 



(28) 



This result is used below in the calculation of the 
Aslamazov-Larkin diagrams, pad may be used to obtain 
the heat current in the TDGLLL3 (again, note the appear- 
ance of the Cooper pair charge, — 2e). To be precise, 
Eq. (28) is the heat current associated with a pairing 
field configuration A = e ' iq ' r ~' ,wt , using [cf. Eq. (|2T 



J Q = -T)V 



<9A* 
~dt 



VA- 



(29) 



Previously, the heat current vertex 3 Q was considered 
by several authors, beginning with the work of Caroli and 
Maki.Ea However, Eqs. ([l9]) and ( p8| ) do not appear la 
have been obtained previously with the correct factor.E3 
The same result may of course be obtained by an explicit 
(but morCpCumbersome) calculation of the heat current 
vertex J Q E3 However, in addition to their being more 
straightforward, the arguments presented in this section 
have the advantage of being very general, applicable to 
arbitrary disorder strength and range. [Note that the 
important ingredient used to obtain Eqs. (|24|), (|25|), and 
( pf ) is just the absence of inelastic scattering.] Finally, 
we have shown here that Eq. ( |l9| ) holds to linear order in 
q; we note that the argument may be extended to higher 
orders in q as well. 



where the function / results from integration over all 
internal momenta in the diagrams. The important point 
here is that the function / that appears in Eq. (|24|) is 
identical to the one that appears in Eq. (p6|). Eq. (]26|) 
may be rewritten as 



J Q = qj^ \ic m ;pj /(e m , -e m + uj m ) 

+ q^^r/(fm,-£m + w m ). (27) 



Here, the first term on the right hand side can be shown 
to vanish using Eq. (|25|) . On comparing the second term 



V. CALCULATION OF THE 
ASLAMAZOV-LARKIN DIAGRAMS 

In this section we calculate the Aslamazov-Larkin con- 



tribution to i 



The starting point is the expression for 



the Aslamazov-Larkin diagrams depicted in Fig. |4|. To 
leading order in momentum and energy, the current ver- 
tices depend only on momentum and energy flowing from 
one fluctuation propagator to the other, as in Eqs. ( p^ ) 
and (pH). It follows that the Aslamazov-Larkin contri- 
bution to the current correlator A [see Eq. (||)] is given 
by 
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A-(Q,zfi m ) = / ^3 ( 3 °) 

x | J^fe + Q) Jy(q y ) Jy{q y ,iu m + i(q, iu; m ) L(q + Qx, iw m ) L(q + Qx, 

+ J^(g x ) Jy(q y , iu m + i(q, iw m ) i(q, iw m + iO m ) i(q + Qx, iw m + iO m )|. 

Following the standard procedure, the sum over Matsubara frequencies may be expressed as an integral over the 
contour presented in Fig. ^ The resulting expression, after the analytic continuation zf2 m — * f2 + iO, is given by 



x je(q x + Q)j;(q v )\jQ(q y ,u + n/2)L R (q + QSt,u J + n)Im i R (q, u)L R (q + Qx, u) 



+ J®(q y ,tO - n/2) L A (q,w - II) i A (q + Qx, w - fi) Im L H (q + Qx 



+ ^(<fr) <^ w + 0/2) L*(q, w + n) L R (q + Qx,w + fi) Im L*(q, w) 



(31) 



J Q ( 9y , w - fi/2) L A (q, w - 0) Im L R (q, w)L H (q + Qx, u) 



Here, n(w) = | coth(w/2T), and L R (q, u) and L A (q, w) 
are the analytic continuation of L(q, ioj m ) on the two 
sides of the cut at Imw = 0. The main contribution to 
the integrals is from small wavevectors and frequencies; 
to leading order in T — T c , n(u>) « T/u>, and the fluctu- 
ation propagator L and current vertices J e and J Q are 
given by Eqs. (|l3|), (p2|), and ( |28| ) respectively. Next, the 
expression is expanded to linear order in Q and f2, the 
integrals are calculated, and using Eq. (|5|), we obtain for 
two and three dimensions, 



Jl 
E 



e 2 TB r] 
2ttc e 



e 2 TB 



for 2D, 



for 3D. 



(32) 



With the identification of the superconducting coherence 
length £(T) = this result is identical to the one 

obtained by considering the Gaussian fluctuations in a 
stochastic TDGlMeI 



As discussed in Sec. II A, it is necessary to subtract the 



magnetization current j^ ag = cM x E from this result to 
obtain the correct transport response. The correspond- 
ing contribution of superconducting fluctuations to the 
magnetization is given by 



M = 



e 2 TB Tj 
37rc 2 e 

e 2 TB 

67TC 2 



for 2D, 



for 3D. 



(33) 




q + Qx, uj m + Q„ 



q, w„ 




q + Qx, lj. 



FIG. 4: Aslamazov-Larkin diagrams contributing to jy / EB 
[See Eqs. (|^)-(^)]. The wavy lines correspond to the fluc- 
tuation propagator L; electric and heat current vertices are 
indicated in the figure. 



We note that in the Aslamazov-Larkin calculation the 
magnetization currents contribute two thirds of the total 
current in both two and three dimensions. 

The final result for the Aslamazov-Larkin contribution 
to a xy is obtained after the subtraction of the magneti- 
zation currents. The result is given in Eq. (0), where we 
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FIG. 5: Contour in the complex uj plane used for expressing 
the sum over Matsubara frequencies uj m in Eq. (BOI) as an 
integral, leading to Eq. (|il|). The contour runs along the cuts 
of L((\,uj) and -L(q, u> + ifl m ), and is closed at a distance from 
the origin which is taken to infinity. The dots are the poles 
of n{ui) at iui m . 

introduce back Ti, and present the result in terms of the 
coherence length, £(T) = \Jrj/e. 

VI. MAKI-THOMPSON AND DENSITY OF 
STATES TERMS 

In this section we consider the Maki-Thompson and 
density of states diagrams [diagrams (b)-(m) in Fig. |). 
We show here that these terms are less divergent than 
the Aslamazov-Larkin diagrams as T — > T c . In addition 
a similar conclusion can be drawn for the magnetization. 
We thus find that in the microscopic calculation of a xy 
the Aslamazov-Larkin result, Eq. (||), is the most diver- 
gent contribution as T — > T c . 

We note that a similar result may or may not hold for 
different transport coefficients, and each case must be 
examined separately. Indeed, in the case of the conduc- 
tivity in three dimensions the Maki-Thompson contribu- 
tion has the .-same divergence as the Aslamazov-Larkin 
contribution^ In two dimensions the Maki-Thompson 
diagram diverges at any temperature in a naive calcu- 
lation, a divergence which is regularised by introduc- 
ing a pair breaking mechanism. ESoEj After this regu- 
larization, the Maki-Thompson contribution has a form 
which is somewhat different than the power law of the 
Aslamazov-Larkin term, and its ultimate divergence is 
only logarithmic. Nevertheless, this is an important mi- 
croscopic contribution to the conductivity (except very 
close to T c ). In contrast, a different situation may hold 
for other transport properties. For example, Niven and 
Smith have recently shown that the contribution of the 
Maki-Thompson and density of states teraas for the ther- 
mal conductivity does not diverge at T c .Eil 

The method that we use here to obtain our result is 
that of power counting, applied to each of the diagrams 
independently. We thus avoid the explicit calculation 
of the diagrams, which would be needed if sub- leading 



terms are desired. In each of these diagrams, after each 
of the electronic blocks is calculated, the structure that 
remains is that of an integral over momentum and energy, 
with the integrand being composed of fluctuation prop- 
agators and electronic blocks. We apply power counting 
arguments to this integral to find the dependence of each 
diagram on the reduced temperature e. (This procedure 
does not exclude the possibility that the coefficient of this 
power is actually zero and that the power of the diagram 
is therefore lower, nor does it exclude the possibility that 
the diagram is identically zero.) 

For the purpose of clarification, we begin by consider- 
ing the Aslamazov-Larkin diagram [Fig. |^(a)] , which was 
calculated explicitly in Sec. M. The power counting is 
thus applied to Eq. ( |3l| ) ; we now count powers of e in the 
integral explicitly. In the integrand, there are three fluc- 
tuation propagators L (contributing a power e _1 each), 
two electric current vertices J e (e 1 / 2 each), and one heat 
current vertex J Q (e 3/2 ). To obtain j^/EB [see Eq. (§)], 
the integral is expanded in external frequency f2 (e _1 ) 
and external wavevector Q (e -1 / 2 ). The integration over 
momentum gives another e d ^ 2 , while there is no contribu- 
tion associated with the integration over energy [because 
of the n(uj) factor]. Accounting for all contributions, we 
obtain a divergence of e d / 2 ~ 2 . Similar arguments give 
an identical result for the divergence of the magnetiza- 
tion, giving a^y oc e rf / 2 ~ 2 ; in agreement with our exact 
calculation, Eq. (||). 

We consider next diagrams (b)-(g) in Fig. ||. In these 
diagrams, the number of fluctuation propagators is one 
less than in the Aslamazov-Larkin diagram. If any of 
these diagrams is to be as divergent as the Aslamazov- 
Larkin diagram, then this loss of a power of (T 1 must be 
compensated (as it is in the case of the Maki-Thompson 
diagram in the conductivity). A power of e _1 is re- 
gained when considering the electronic block of the di- 
agram which has two microscopic current vertices in it, 
provided they are in the same direction. Here, this will 
only occur for diagrams (d) and (e), in which the two 
vertices which are in the same block are if and j®. In- 
deed, in the Aslamazov-Larkin diagram, the vertices Jy 
and J® contribute e 2 . In diagrams (d) and (e), the block 
containing the vertices jy and contributes only one 
power ofe.ll Moreover, for these two diagrams to diverge 
as the Aslamazov-Larkin diagram, the expansion in ex- 
ternal wavevector Q and external frequency Q should give 
powers of (T 1 !" 1 and respectively, as it does for the 
Aslamazov-Larkin diagram. The expansion in indeed 
gives a power of e _1 , since the external frequency appears 
explicitly in the fluctuation propagator (and also due to 
a diffusive pole as in the case of the Maki-Thompson con- 
ductivity diagram). However, the important point in this 
analysis is that in diagrams (d) and (e) the expansion in 
external momentum Q does not gain a power of e -1 / 2 , 
but instead accounts for e 1 / 2 in the power counting. The 
reason for this is that the external momentum, which 
flows from the jy vertex to the vertex, does not flow 
through the fluctuation propagators of the diagram. The 
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expansion in external momentum is thus limited to the 
electron block which includes these two vertices, where 
it is straightforward to check that expansion in the ex- 
ternal momentum leads to a power of e 1 / 2 . Finally, di- 
agrams (h)-(m) involve only one fluctuation propagator 
and will clearly be less divergent than the Aslamazov- 
Larkin term. 

We have demonstrated by power counting arguments 
that while the Aslamazov-Larkin diagram, Fig. ||(a) , di- 
verges as e d / 2_2 ; all other diagrams in Fig. g are less 
divergent as T — > T c . Similar arguments hold for the 
fluctuation contribution to the magnetization, and hence 



the critical regime (which for low-temperature supercon- 
ductors is very narrow, as expressed by the Ginzburg cri- 
terioncl) . Additional microscopic terms become irrelevant 
in this regime. 

On the other hand, further away from T c in the re- 
gion where the microscopic calculation is valid, additional 
microscopic contributions may arise (as they do for the 
conductivity). To reiterate, for the transverse thermo- 
electric response we find that they are less divergent then 
the Aslamazov-Larkin contribution as T — * T c . Not in- 
vestigated here is the possibility that these normal state 
corrections vanish in the case of particle-hole symmetry, 



for a xy . In view of this conclusion we do not calculate as does the Drude result for a xy (to emphasize, a result 



the Maki-Thompson and density of diagrams explicitly. 
[We note that diagrams (b)-(g) in Fig. [2| as well as the 
Maki-Thompson and density of states diagrams in the 
calculation of the magnetization, will have a logarithmic 
divergence in two dimensions as T — > T c .] 



VII. CONCLUSIONS 

The main result of this paper is that the leading contri- 
bution in the microscopic calculation of a xy arises from 
the Aslamazov-Larkin diagrams, which correspond to the 
contribution of Gaussian fluctuations in the stochastic 
TDGL. The Maki-Thompson and density of state terms 
are less divergent as T — > T c . In concluding this paper, 
we comment on several aspects of this result. 

It is well known that in calculating the contribution 
of superconducting fluctuations to the conductivity, the 
Aslamazov-Larkin contribution corresponds to the Gaus- 
sian contribution of a stochastic TDGL.EJtl In establishing 
the correct form for the heat current vertex in Sec. IV, 
we verify this correspondence also for thermal transport. 

There are two directions in which to extend the cal- 
culation beyond this approximation: by considering the 
additional microscopic contributions (as we did in this 
paper), or by going beyond the Gaussian approximation 
in the stochastic TDGL (cf. Ref. |f). We would like to 
emphasize that these approaches are of a very different 
nature, and their regime of validity is also different. 

The stochastic TDGL is traditionally understood as 
the model for the critical dynamics of a superconduc- 
tor (model A in the classification of Hohcnberg and 
HalpcrmEj). As such, the TDGL should give the rele- 
vant contribution as the temperature approaches T c in 



not required by symmetry). 

To connect the microscopic approach with the criti- 
cal dynamics, one may expect that in the microscopic 
approximation as the temperature is lowered, the Maki- 
Thompson and density of states terms would become 
less important as the behavior becomes governed by the 
stochastic TDGL. Fox-u. xy this clearly occurs in the mi- 
croscopic calculation!!-] In conclusion, this work provides 
further justification for using the TDGL also as the tem- 
perature increases away from T c into the Gaussian regime 
(which is the approach taken in Ref. 

Interest in the Nernst signal has grown recently due 
to the measurements in high-temperature superconduc- 
tors, whepe,|the fluctuations signal can be observed well 
above T c Jl!3t!J On the other hand, the contribution of 
superconducting fluctuations to the Nernst signal is yet 
to be observed in a low-temperature superconductor, for 
which the BCS microscopies considered in this paper are 
applicable. As discussed in the Introduction, we expect 
the fluctuation tail to be observable in the Nernst signal 
of a suitably chosen superconducting thin film. It would 
certainly be of interest to verify this prediction experi- 
mentally. 
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